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Abstract —This paper proposes a novel approach to phase- 
noise compensation. The basic idea is to approximate the phase- 
noise statistics by a finite number of realizations, i.e., a phase- 
noise codebook. The receiver then uses an augmented received 
signal model, where the codebook index is estimated along with 
other parameters. The realization of the basic idea depends on 
the details of the air interface, the phase-noise statistics, the 
propagation scenario and the computational constraints. In this 
paper, we will focus on a MQAM-OFDM system with pilot sub¬ 
carriers within each OFDM symbol. The channel is frequency 
selective, fading and unknown. A decision-feedback method 
is employed to further enhance performance of the system. 
Simulation results are shown for uncoded and coded systems 
to illustrate the performance of the algorithm, which is also 
compared with previously employed methods. Our simulations 
show that for a 16-QAM coded OFDM system over a frequency 
selective Rayleigh fading channel affected by phase noise with 
root-mean-square (RMS) of 14.4 degrees per OFDM symbol, the 
proposed algorithm is 1.5dB from the ideal phase-noise free case 
at a BER of 10 -4 . The performance of the best reference scheme 
is 2.5dB from the ideal case at BER of 10 -4 . The proposed scheme 
is also computationally attractive. 

Index Terms —OFDM, phase-noise, BER and SNRs, codebook, 
fading channel, channel estimation, pilot sub-carrier. 


I. Introduction 

Orthogonal frequency division multiplexing (OFDM) is 
adopted by most of the current and future telecommunica¬ 
tion standards for high-rate data transmission, particularly in 
wireless communication systems. Its resilience to multipath 
channel fading, the spectral efficiency it provides as well as the 
simplicity of the equalization, has enabled OFDM to remain 
the most popular modulation scheme. However, OFDM is 
known to be sensitive to various hardware imperfections, the 
Dirty-RF effect Ql , originating in the transceiver hardware. 

A. Previous Works 

The phase-noise (PHN) phenomenon and its underlying 
effects on various OFDM systems were studied extensively 
in iUl. The effect of PHN on OFDM systems is classified 
into two components. The common phase error (CPE), which 
rotates all the sub-carriers in one OFDM symbol by a common 
phase distortion, and the inter-carrier interference (ICI), which 
arises due to the loss of orthogonality between each sub¬ 
carrier. Frequency domain approaches to PHN estimation and 
compensation mainly deal with the CPE and ICI components 


separately, while time domain approaches attempt to com¬ 
pensate for both jointly. It has been shown that a significant 
improvement in performance can already be achieved with 
CPE correction only by treating the ICI term as an additive 
Gaussian noise [sj-11]. However, CPE correction only is not 
always sufficient for high rate transmission, therefore ICI 
compensation is necessary. 

Most ICI compensation techniques employ decision- 
directed feedback (DD-FB) for frequency domain PHN esti¬ 
mation Q 12-H. A vast majority of works assume a known 
channel frequency response iH TtH [l5l - fl9ll . For unknown 
channel, joint channel and PHN estimation is considered in 
113] where an ICI reduction scheme over a Rayleigh fading 
channel is presented in which the PHN process within an 
OFDM symbol is modeled as a power series. Although the 
method presented showed a significant bit error rate (BER) 
improvement, the cost in computation is large. A less complex 
method is presented in |14] where the estimate of channel and 
CPE at pilot subcarriers are interpolated to obtain channel 
frequency response and CPE followed by a DD-FB loop to 
estimate ICI components. 

A non-iterative compensation scheme in which the CPE 
between consecutive OFDM symbols are interpolated linearly 
to estimate the time varying PHN is presented in E3 0. 
In [20], a method of suppressing the ICI is presented by 
linearly combining the cyclic-prefix and the corresponding 
OFDM samples. The linear coefficients are obtained such that 
the ICI power is minimized. Other methods using the Bayesian 
framework for joint channel and PHN estimation have been 
proposed in 121 42411 . Soft-input maximum a posteriori and 
extended Kalman smoother are proposed in 121]. In [22- 


2411 . Monte-Carlo methods are employed to approximate the 


posterior probability distribution of the unknown quantities 
for PHN tracking and channel estimation. Although they 
are shown to provide good BER performance, Monte-Carlo 
methods, such as particle filters, are known for their numerical 
complexity. The main interest in the paper lies on PHN 
compensation over unknown channels using pilot sub-carriers. 
However, the proposed method also uses DD-FB loop for 
improved performance. For the known channel the proposed 
method is benchmarked against @]. For known and unknown 
frequency selective Rayleigh fading channels the methods in 
ni3L [l4] are also used as benchmarks. The problem setting in 
which the methods presented in 0, dill and are employed 
is identical to the setting the proposed algorithm. 
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B. Contributions 

In this paper, we employ a simple and novel scheme for 
PHN cancellation which is shown to have superior perfor¬ 
mance in terms of BER compared to previously proposed 
methods. The proposed technique can be applied to any 
assumed distribution of the PHN process, 9{n). Additionally, 
no approximation based on the magnitude of the PHN is taken 
into account, e.g. e~ 1 + j sin(0(n)), as is common in 
some previous works. 

A codebook of K vectors representing a set of trajectories 
which aim to closely match with the PHN realization is used. 
An uncountable set of possible PHN realizations is represented 
by a codebook with K quantized trajectories, which are stored 
at the receiver^ The complexity of the algorithm is therefore 
dependent on the number of trajectories in the codebook, K. 
For a moderately small sized codebook, it is shown that the 
impact of the PHN is significantly reduced (e.g. with K = 27, 
there is 60% reduction in terms of the effective PHN mean 
square error (MSE), see Section Ull-Ck 

In this paper, we focus on a solution where the tra¬ 
jectory that minimizes the Euclidean distance between the 
constellation of the received symbols at the pilot positions 
and the constellation of the known pilot symbols is chosen. 
Additionally, a DD-FB technique is employed such that both 
the estimate of data symbols from the channel decoder and 
pilot symbols are used to compute the Euclidean distance. The 
proposed technique can also be used for combined channel 
estimation and PHN compensation schemes, in which MMSE 
based channel estimation is employed for each of the vectors 
in the codebook. The channel estimate corresponding to the 
trajectory in the codebook which best approximates the PHN 
realization is chosen. We will show by the simulation that, 
for AWGN and fading channels with a known frequency 
response, the proposed scheme outperforms previously pro¬ 
posed PHN compensation techniques. In this scenario, the 
proposed method does not employ DD-FB while the reference 
schemes require previously detected symbols for ICI estima¬ 
tion. Simulation results are also shown for unknown channel 
frequency response, in which the proposed method employs a 
DD-FB technique for combined channel estimation and PHN 
compensation. It is shown that the method presented provides 
improved performance compared to previously proposed meth¬ 
ods. Additional results based on a PHN process modelled as 
Ornstein-Steinbeck is also shown to illustrate the applicability 
of the proposed method for a PHN model other than Weiner 
process. 

The paper is organized as follows. Sec. [TT] introduces the 
system model and parameters while Sec.QII]presents design of 
the codebook for Weiner PHN along with some performance 
analysis. In Sec. [IV] details of implementation of the code¬ 
book for PHN compensation together with channel estimation. 
Simulation results are shown in Sec. [V]and a computational 
analysis of the proposed algorithm as well as some of the 

'The technique is similar to the design of vector quantization codebook 
where a set of n vectors from some m-dimensional space is efficiently 
represented by a codebook with a smaller set of vectors from m-dimensional 
space j251. 


reference methods is given in Sec. [VI] We end with concluding 
remarks in IVIII 


II. System Model and Proposed Approach 


We consider an OFDM transmission system over a fading 
time varying channel with L multipath taps employing QAM 
modulation technique. The OFDM system is assumed to have 
FFT size of N of which N p are pilot subcarriers and with a 
cyclic prefix of length N cp such that N cp > L. To obtain the 
m-th OFDM symbol of duration T, a stream of data bits is 
divided to N—N p groups of M- bits which are mapped to 2 M - 
QAM complex symbols, S m (k). Assuming perfect frequency 
and timing synchronization, the baseband representation of the 
received time domain signal, y m (ri), is given by 

y m {n) = e j9mM f ^2 hm{Z)s m (n - ?)) + w m {n) (1) 


i -o 


where s m ( n ) is the N point inverse FFT (IFFT) of the complex 
transmitted symbol, S m (k), h m (i) is the fading complex 
channel pulse response with L propagation paths, w m (n ) is 
a zero mean complex circular Gaussian channel noise with 
variance cr^; and 9 m (n) is the PHN sample at the receiver at 
time index n of the m-th OFDM symbol. The received signal 
is sampled at a frequency f s = (N + N cp )/T, where T is the 
OFDM symbol duration. 

After removing the cyclic prefix, an FFT operation is 
applied to the received signal such that the fc-th subcarrier 
is given by 


N-1L-1 

Y m (k) = -= Y, E h m(0Sm(n - i) e ^{n)-2,kn) 
v N n —0 £=0 

n =0 

N-l 

= S m (k)H m (k ) A m ( 0) + ^2 Sm(l)H m (l)A m (l - k) 

CPE \2k 

' -V-' 

ICI 


+ w m (k) 


( 2 ) 


where A m (k) and W m (k) are respectively the FFT of eJ 9m< ' n ' > 
and that of the additive complex Gaussian channel noise w[n). 
H m (k) is the instantaneous frequency response of the multi- 
path channel which is assumed to be constant within one 
OFDM symbol. Figure [j] shows the block diagram of the 
system model. Equation © can further be written in matrix 
form as 

Y m = A m S^H m + W m (3) 

where = diag(S , m (0),.., S m {N — 1)) is an N x N matrix 
representing diagonalization of the symbol vector S m ; and the 
N x N matrix A m is a unitary circulant matrix containing 
DFT coefficients of e^ 9rn ^ such that A m (k, l) = A m (l — k). 
H m and W m are vectors of the channel frequency response 
and the complex Gaussian noise respectively. Estimation of the 
realization of the PHN process, 9 m = [0 m (O), • • ■ 9 m (N— 1)] T , 
would enable compensation of CPE as well as ICI. 
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Figure 1: Block diagram of the OFDM system. 


A. Approximate Model 


the criterion 0 as 


We propose a method where the PHN realization at the m-th 
OFDM symbol is approximated by one of a set of K vectors, 
where = [4> k {0),<j> k (l), -■ ■ , <j> k {N-l)] T plus an 
initial phase offset, t/> m € [0, 2n). Associated with each vector, 
<j> k , is a corresponding unitary circulant matrix containing DFT 
coefficients of {e- 7 ^ fe ^}^J 0 1 denoted by A. The signal model 
in 0 is then modified as 

Y m = A km S^H m + W m (4) 

which is an approximation of ©, where the PHN DFT 
coefficient matrix, A m , is now replaced by e^ m A km , which 
is a function of an index, k m , that takes on a value between 1 
and K with equal probability. The framework can be extended 
to unequal probabilities but this is left for future work. The 
random variable ip m , which is assumed independent of k m , 
is uniformly distributed and models the initial phase of the 
symbol. Both k m and i/j m are independent from symbol to 
symbol. The design of the codebook {(f> k } k=1 is yet to be 
discussed in a later section. It should be noted that this 
approximate model is only used for the purpose of deriving the 
algorithm. Data from the original model in (J3]» is used when 
testing and validating the proposed algorithm by simulation. 


B. Cost Function 

Any estimator derived in the absence of PHN can be written 
as 

H = argminC(Y 1; .. ., Y g ), (5) 


fl = arg 


Q,ki 


min C(e 


—3i >i A —1 


A^Yr, 




q A 


( 6 ) 


for the model (0. This includes, for instance, maximum- 
likelihood estimators where now the estimator in © is also 
the maximum likelihood estimator for the model 0}. An 
example of a use of the proposed new estimator would be 
in an OFDM system with an equalizer and Viterbi decoder 
and a pilot subcarrier. In such a system, the optimization 
function could be the sum path metric. In this case, 0 implies 
that the entire receiver is run for all possible hypotheses of 
ki... k q (of which there are K q ), and selects the one with the 
lowest path metric (in each symbol optimizing over i/j). This 
implementation is generally not realistic. A practical criterion 
is hard decoding of the M-QAM symbols of OFDM symbols, 
where the criterion function separates into a term for each 
symbol individually. In this case, the criterion function is the 
Euclidean distance between estimates of the received signal 
point, S m , at the output of the equalizer and the nearest point 
in the M-QAM constellation. Here, it is realistic to sequence 
through a quite large size K. Selection of the codebooks is a 
design problem still to be considered. 

For the case where the Wiener PHN model is used, a 
systematic codebook design is presented in the next subsec¬ 
tion. For other PHN models different codebook designs may 
perform better. However, the design for Wiener PHN may 
work well for other PHN models as well. This is illustrated 
in Section IV-B5I where results based on a PHN modelled as 
Ornstein-Uhlenbeck process is presented. 


where O is the estimated parameter, C is some form of 
optimization criterion and q is the number of OFDM symbols. 
It should be noted that A km and its inverse A ; ( 1 are unitary 
circulant matrices. As a result. At is invertible and AT 1 W,„ 

m 0-771 

is statistically identical to W m . We may therefore generalize 


III. CODEBOOK DESIGN FOR WEINER PHN 

The model introduced in jH] and commonly employed in 
literature, where the PHN process is described by a Wiener 
process, is used. Under the Weiner model, the PHN sample at 
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Figure 2: A single Wiener PHN realization (blue) with crj = 2 7r (0.01/TV) 
where N = 64 and its approximation according to © for J = 8. 



Figure 3: Dividing the Gaussian sample space into Q regions of equal areas, 
where Q is an odd integer. 


the m-th OFDM symbol is given by 

mN -\-n n 

0 m {n) = ^ e(i) = V> TO + y^e(Q (7) 

i —0 1=0 

for n = 0, 1,. .., N — 1, where e(Z) is a Gaussian random 
variable with zero mean and variance cr 2 = 2nf3T/N (in 
radians 2 ), and f3T denotes the rate at which the PHN variance 
grows in one OFDM symbol. The objective is to construct a 
codebook containing a finite set of vectors which approximate 
the trajectories of possible Wiener PHN realizations within 
one OFDM symbol. The subscript m in 6 m {n) and ip m [n) is 
dropped for notational simplicity. Let the N samples of the 
OFDM symbol be divided into J segments with L = N/J 
samples per segment. Then 9{n) is defined as a process which 
is constant in each segment. The value of the constants, A j, 
in each segment are set as the sample average of 9(n) in 
the corresponding segment. The PHN process is therefore 
approximated by 9(n) given by 

1 i_1 

§(n + jL) = \ j = I '520(t + jL) 

. f-'U+r 

= ^ E E eW+V’ 

f =0 z —0 


for n = 0, • ■ • , L— 1 and j = 0, • ■ • , J—1. This is illustrated in 
Figure[2]for <r 2 = 2ir(0.01/N), N = 64 and J = 8, where the 
samples within each segment are approximated by the sample 
mean of the elements in the segment. The process, 9(n), is 
similar to a Weiner process where each step, Xj = A / + i — A j, 
occurring every L samples, is given by 

Xj = I E i 0 ^ + + 0 - + *)) 

, L -1 l+(j+l)L ^ 

= i E E «(0 

^=0 i=e+jL +1 

where the second equality follows from ©. The increments, 
Xj, are identically distributed Gaussian random variables with 
zero mean and variance a 2 r given by 

, L — l £+(j+l)L L — l v+(j+l)L 

al = E[X*] = -J2 E E E E H0e(k)} 

£=0 i=£+jL -\-1 v=0 k=v+jL -\-1 
_ (2£ 2 + 1) 2 
3 L CTe 

( 10 ) 

where the auto-correlation E[e(£)e(k)\ = <r 2 S(k — £) as well as 
arithmetic and geometric series are used to obtain the second 
equality. Note that for L = 1, that is when the number of 
segments is equal to the number of samples, J = N, then 
(j 2 = of. The approximated Weiner PHN process can be re¬ 
written in terms of Xj as 

i 

9(jL + n) = Aj_i + Xj + ip = ) ' Xi ij) (11) 

e=o 

where A_i = 0 and Xq = 0 for j = 0. 

Since the Gaussian sampling space is unbounded, the ran¬ 
dom variable X is instead represented by a set of quantized 
samples, {xi}f =1 , which represent Q regions, {7?.,;}^, in 
the Gaussian pdf with zero mean and variance <j 2 . This 
implies that the increments at the consecutive segments, 
j = 1, • • ■ , J — 1, are defined by a set which define 

Q possible trajectories. That is, the second segment, j = 2, 
has Q possible increments which set Q possible trajectories. 
In the next segment, there will be Q possible increments on 
each trajectory from the previous segment such that there are 
Q :i ~ 1 trajectories at the j-th segment. Therefore, at the J-th 
segment, there will be K = Q J ~ 1 trajectories representing 
quantized paths of the random walk process each of which 
are set as a vector entry in the codebook, 4>^ =1 . The number 
of code vectors, K, is therefore determined by the predefined 
set of quantization regions, Q, and the number of segments 
the symbol is divided into, i.e. J. 

A. Defining Quantization Regions 

The quantization regions can be defined as a set of regions, 
{Ki}? =1 , which divides the Gaussian sample space into Q 
regions. For ease of development, the sample space is divided 
into a set of regions with equal areas (i.e. equiprobable). The 
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Figure 4: A single Wiener PHN realization (blue) with a^ = 2tt(0.01 /N) 
where N = 64 and a finite set of K = 2187 trajectories,* for J = 8 and 
Q = 3, which corresponds to possible trajectories that the PHN realization 
might take for equiprobable quantization. 


regions are defined as 


P 


X G Pi 


= ■•• = P 


X G P, 


Q 


Q 


( 12 ) 


where P[A” G Pi] implies the probability that the random 
variable X is in the region P,. For the case when Q is odd, 
the sample space is partitioned by Q + 1 data points denoted 
by {—oo, X- r ,... , X-i, xi, . .., x r , oo}, as shown in Figure 
0] where X-t = — xi, t = 1,... ,r and r = (Q — l)/2. The 
regions at the tails of the Gaussian pdf are given by TZi = 
[—oo, x_ r [ and Pq = [x r , oo[, while the regions in the middle 
are bounded by [xt,xt+i[. The bounding points, xg are given 
by 





0-1 

2 


03) 


where erf _1 (-) is the inverse error function. On the other hand, 
when 0 is even, the sample space is partitioned by 2r +1 data 
points denoted by {—oo,X- r ,- ■ ■ ,X-i,xo,Xi,- ■ ■ ,x r ,oo} 
where xq = 0 and X-t = —Xu for t = 0, ••• ,r where 
r = Q/2. These points are given by 



The plot in Figure 0] is shown as an example to illustrate the 
process for a set of equiprobable quantization regions where 
0 = 3 and J = 8 while Figure 0 shows the trajectory that is 
closest to the PHN realization. The regions can also be defined 
as a set which divides the Gaussian sample space uniformly. 
This implies that the codebook trajectory would have unequal 
probability. The development of this case is left for future 
work. 


B. Quantization Points 

Given the quantization regions, P = {Pi,...,Pq}, the 
quantization points, {xi,... ,xq}, which represents each re¬ 
gion is defined as the mean point within the region. That is. 


Xi = E 


X\X G Ri 


S Ri xfx{x)dx 
Jr. fx(x)dx 


(15) 



Figure 5: The trajectory that best matches the PHN realization. 


for i = 1 ,... ,Q, where E[-] is the expectation operator and 
f(x) is the Gaussian distribution function. It can be shown that 
( fT5l > provides the set which results in a minimum squared error 
quantization of the random variable X with a pdf f(x), (25]. 
Moreover, the expected value of quantization error q er ror = 
x — x is zero, i.e. E[x] = E[x]. For the Gaussian pdf with 
zeros mean and variance cr^ = (2 L 2 + l)/3 LaJ ( IT5] > is given 
by 


f = -o-x / exp(-0.5xf +1 /crg) - exp(-0.5xf/ag) \ 

O.bV^TT V erf(xi+i/\/2cr x ) - erf(xj/\/2cr x ) J 

(16) 

It can be seen that the codebook will contain an all zeros entry 
when Q is odd. This is attractive, since it means no correction 
will be done if the receiver would have no phase noise and the 
SNR is high. The plot in Figure 0] employs (03) for Q = 3 to 
define the quantization regions and 06) to obtain the quantized 
point for each region. 


C. MSE Analysis 

In this section, the MSE for the designed codebook is 
analysed. The objective is to provide a good indication of 
the performance of the designed codebook for a given J and 
Q. The MSE between a random PHN realization and the 
approximate model introduced is given by 

MSE = E 0(n) [miny 1 9(n) - ip- <p k (n) | 2 ] (17) 

k,ip z ' 


which is difficult to evaluate analytically. Therefore the MSE is 
evaluated by simulation. However, an approximate expression 
for the MSE, derived in Appendix |A] is given by 


MSE = 


(N+J)(N-J) 2 


6 J 


o e + L(J - 1 )a 


(18) 


with a q given by 06). 

Table U shows the MSE values obtained by Monte-Carlo 
simulations (denoted MSE s ) normalized by the MSE for the 
CPE correction only (i.e. MSE for J = 1). The number of 
PHN realizations employed were 5000 for every codebook 
with various combination of J and Q resulting in K = 
Q J ~ 1 . The approximated theoretical MSEa values computed 
according to 08) are also given. It can be seen that the 
analytical and simulation results are similar indicating that 
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Figure 6: Block diagram of the coded OFDM receiver with PHN compensation and channel estimation. 


Table I: Approximated MSE for N = 64 as is given by GU and simulated MSE denoted by MSE s for various values of Q and J. The values are normalized 
by the MSE value for CPE correction only given by MSE = (N — 1 )(N + 1)(T^ /6. 



K 

<2 = 2 
MSE s 

MSE a 

K 

Q = 3 
MSEs 

MSE a 

K 

<2 = 4 
MSEs 

MSEa 

K 

<2 = 5 
MSEs 

MSEa 

K 

<2 = 6 
MSEs 

MSEa 

J = 

1 

1 

0.9967 

1 

1 

0.9997 

1 

1 

0.9998 

1 

1 

1.0021 

1 

1 

1.0043 

1 

J = 

2 

1 

0.7353 

0.7387 

2 

0.6095 

0.6247 

3 

0.5674 

0.5765 

4 

0.5432 

0.5519 

5 

0.5302 

0.5376 

J = 

4 

8 

0.4479 

0.4309 

27 

0.3488 

0.3446 

64 

0.3076 

0.3082 

125 

0.2878 

0.2895 

216 

0.2770 

0.2785 

J = 

5 

16 

0.4839 

0.3607 

81 

0.2931 

0.2855 

256 

0.2553 

0.2538 

625 

0.2370 

0.2375 

1296 

0.2274 

0.2279 

J = 

8 

128 

0.2631 

0.2330 

2187 

0.1905 

0.1813 

16384 

0.1601 

0.1595 

78125 

0.1487 

0.1482 

279936 

0.1441 

0.1416 


the reasoning used in the analytical derivation in Appendix 
IaI is correct. The MSE analysis assumes that the optimum 
code-book entry is chosen. In order to account also for the 
effect of noise and decision feedback errors, simulations of the 
full receiver for various codebook sizes are performed in Sec. 
m These results confirm that the performance improves with 
increasing codebook size. The choice of codebook size will 
ultimately depend on the computational cost and performance 
requirements for the application at hand. However, it is noted 
that complexity grows rapidly with Q and J, thus moderate 
numbers such as e.g. Q = 2, J = 5 or Q = 3 and J = 4 
seems to provide a good compromise. 


IV. Implementation on OFDM System 
A. Known Channel Response 

The received signal vector after removing cyclic prefix, 
y rn = [llm (0), * * * ,y m {N - 1)] T , is multiplied by each of 
the K trajectories e~^ k such that a new set of K de-rotated 
OFDM signals is obtained at the receiver 

y k m = [e~ jM0) ym( 0), • • • , e-^ N -Vy m (N - 1)] T . (19) 


An FFT operation on each de-rotated signal vector provides 

( 20 ) 


Y fc _ A-iy 
1 m ~ 1 n 


= A- 1 (A m SlH m 


-W m ), 


where the Nx N circulant matrix A ; 7 ' is a matrix of DFT co¬ 
efficients of e^ Jr!>k<J 'K Frequency domain channel equalization 
is applied on Y^ n , and the criterion by which the trajectory 
that best approximates the PHN realization is given by 


p -1 


k* = argmin , 

Ic • ^ 


i =0 


s^ii) - sm 


( 21 ) 


where £ (To, ■ • ■, £p-i) denotes pilot subcarriers, 

is an element of the estimated symbol vector for the fc-th 

trajectory given by 

= Vk 'diagCH™)- 1 ^ (22) 


with complex term r/k that corrects for the effective CPE which 
includes the offset ip m at each symbol, the DC level of the 
codebook vector, <p k , and the accumulated quantization error. 
It is given by 


V k 


Ef = ~o 1 S^Y^/H^) 

EfJo\s m m 2 


(23) 


In a DD-FB loop in which the output of the channel decoder 
can be exploited, non-pilot symbols can also be used such that 
the criterion in ( l23l ) takes into account the decoded symbols 
as well 


N-l 

k* = argmin N 

k ' ^ 

1=0 


S k m (l)-S m (l) 


(24) 


for given by 


Vk 


^~\N —1 
2^i =o 


S* m (l)Y k (l)/H^(l) 


N — 1 
2^i =0 


l^ m (0l 2 


(25) 


The symbol sent to the decoder is then given by the index 
corresponding to Sj'; t . 


B. Combined Channel Estimation and ICI Cancellation 

Since the channel is not always known at the receiver, it 
needs to be estimated. Assuming the channel remains station¬ 
ary within the time period of an OFDM symbol, for a code¬ 
book of size K, there will be K received signal candidates, 
Y^, which are given by ( l20l ) where the channel response, 
H m (k ), is considered to be unknown. The MMSE estimate 
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of the channel frequency response vector, H^, employing the 
/c-th trajectory of the codebook is then given by [ 2611 


= E[ ^[H m (H m ) H ] + ^In) H* 


(26) 


where E[H m (H m ) n ] is the autocorrelation matrix of the 
channel frequency response for the given statistical model of 
the channel and H* rn is the least squares estimate of H k n 
given by 


tt/c _ o-lyfe 

-^■LSjra 1 m 


' Y£( 0) 
5 m (0)’ 


Y£(N~ 1) 


i T 


Sm(N-l) 


(27) 


However, since not all transmitted symbols are known at 
the receiver, N p pilot symbols, fitted evenly among the N 
subcarriers, are employed to obtain a least squares estimate 
of the channel frequency response at the pilot positions. For 
the sake of convenience, the LS estimate at pilot positions are 
denoted by the vector, H^’ ^- Therefore, the MMSE estimate 
of the channel frequency response, for fc-th trajectory in the 
codebook, is initially obtained based on pilot symbols 


H k m = E[H m (HP m ) H } + gl 




-1 


\lk,P 

A LS,m 

(28) 


where (■) p denotes a vector whose elements are positioned 
at pilot subcarriers. A set of K channel frequency response 
estimates, corresponding to each trajectory in the 

codebook is available at the receiver. Equation (l24l> is used to 
determine the trajectory which closely matches the PHN real¬ 
ization, and the corresponding channel estimate is then used to 
demodulate and obtain rough decisions on the data symbols. 
After channel decoding, a decision feedback technique is then 
employed to obtain (l2Tb based on rough decisions on the 
symbols containing pilot and decision feedback symbols, S m , 
which is then used in (l26l > to compute an MMSE estimate of 
Hjjj with improved accuracy. 

For deeply faded channels, estimation accuracy during 
symbol m can be improved by including previously decoded 
symbols, i.e., S m _i, S m _ 2 > ■ ■ ■ , S m _£>+i, in the estimation 
vector. Therefore, taking into account D previously decoded 
OFDM symbols, the MMSE estimator in d28l > becomes 


H k m = E[H m (Hr] £[H(H) 


\Hi 


e k 


IDN H ls (29) 


where H = 

\ffkm TT^m — 1 
L^LS,??! 1 -“^LSjm —1 ’ ' 


[H m , H, 


TT^m — D + l 
’ n LS,m-D+lJ 


Hm—D+l}T and Hls — 

7 On the initial run of the 
m-th OFDM symbol, symbols at pilot subcarriers are used 
in order to compute (l28l t. after which the criterion in (|2IT > is 
used to choose the best match index, k*, for the trajectory 
which closely matches the PHN realization, (f> k , as well as 
the corresponding channel estimate, . Pilot and decision 
feedback symbols can then be used in the next iteration to 
obtain a better estimate, using (l26l > or ( 1291 if using 
L previous OFDM symbols and (l29l > to obtain k*. Joint 
channel equalization and PHN compensation is employed on 
the received signal Y m such that the symbols sent to the 
decoder are given by 


Algorithm 1 The proposed receiver algorithm with no DD-FB 
for uncoded system and known channel 

1: procedure 

2: Received m-th OFDM symbol 

3: for k = 0 : K — 1 do 

4- Y fe — A~ X Y 

x m — -^k 1 m 

5: Compute r/k in ( 123b using pilot symbols 

6: Compute F(k) = |S* (£) - S m (£i)\ 2 

7: end for 

B: Compute k* = argmim- F(k) in ( 1211 to obtain k* 

9: return 

10: end procedure 


oo) 

where k* denotes the index of the most likely trajectory in 
the codebook obtained from (f24t . Figure [6] shows the OFDM 
receiver with the proposed PHN compensation and channel 
estimation scheme. An outline of the proposed algorithm for 
a known channel response without DD-FB loop is presented in 
Algorithm Q] In Algorithm [2] the algorithm for the combined 
channel frequency response estimation and PHN compensation 
scheme is presented for a coded frame of M OFDM symbols. 


Algorithm 2 The proposed algorithm for channel estimation 
and PHN compensation with DD-FB and channel coding 


1 

2 

3 

4 

5 

6 
7 
B 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


procedure 

Received frame 
do Iteration on DD-FB loop 

do for each OFDM symbol in the frame 
for k = 0 : K — 1 do 
Y fc — A -1 V 

x m — 1 m 

if Zeroth iteration then 

Channel estimation using pilots, ( 1281 
Compute r/k in 1231 using pilot symbols 
Detected symbol for each k using 01 
Compute F(k) = l^(fi)- 

Sm(£i )\ 2 

else 

Channel estimation employing ( 1261 
Compute rjk using 01 
Compute F(k) = |5* (&)- 

Sm(£i )\ 2 

end if 
end for 

Compute k* = argminj, F(k) in ( 1211 (zeroth iter¬ 
ation) jor ( 1241 to obtain k* 

Select 

while End of frame 
Decode frame and deinterleave 
Interleave, channel ceding ami modulation frame 
DD-FB frame S = [So, • • • , S m , • • • , Sm-i] 
while End of iteration 
end procedure 


V. Simulation Results 

In this section, simulation results are given to show the 
BER performance of the proposed technique under Wiener 
PHN. Two sets of simulation results are presented where in 
the first set, the performance of the algorithm in Algorithm 
0] for an uncoded system is presented under AWGN channel. 
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Figure 7: BER performance of proposed algorithm compared to 0 for 
16-QAM and 64-QAM uncoded OFDM under AWGN channel and PHN. 
Corresponding performance plots for an ideal case (i.e. AWGN channel no 
PHN) is also given. BT = 0.01. 


In the second set, performance of the algorithm in Algorithm 
0 is shown for a Rayleigh fading channel. Simulation results 
are presented for 16-QAM and 64-QAM modulation schemes 
with N = 64, N cp = 16 and N p = 8 pilot subcarriers. The 
codebook size for the results presented is K = 27, in which 
Q = 3 and J = 4, which from the analysis in Sec. IIII-CI 
is a good compromise between complexity and performance. 
Perfect timing and frequency synchronization is assumed at 
the receiver. 


A. Uncoded System, AWGN Channel 

Considering the AWGN channel, the proposed Algorithm Q] 
was compared with the PHN compensation method presented 
in 11311 . employing 3 iterations in the DD-FB. In the proposed 
method no DD-FB iterations are employed. The rate of growth 
of the PHN variance /3T is set to 0.01, and 16-QAM OFDM 
and 64-QAM OFDM were considered without channel coding. 
The result is given in Figure [7] For the 64-QAM system 
(dashed line), the proposed method is 2.5dB from the ideal 
PHN free case at a BER of 10 -2 , while the method in 1131 
is 9dB from the ideal case, i.e., an improvement of 6.5dB. In 
the 16-QAM system (solid line), the proposed method is 2dB 
from the ideal PHN free case at a BER of 10~ 3 , while the 
method in [[13] is 3dB from the ideal case. 


B. Coded System with Fading Channel 

A 10 tap Rayleigh fading channel with parameters specified 
in Table HU is used. Moreover, a 1/2 rate convolutional code 
with a constraint length of 7 is used on the input bit stream, 
which is then passed on to a bit interleaver over 20 OFDM 
symbols. A Viterbi decoder is used at the receiver. The PHN 
variance growth rate, f3T, is set to 0.01 for 16-QAM OFDM 
such that the RMS in degrees is 180y/2nj3T / n = 14.4°; while 
for 64-QAM, f3T is set to 0.005. Higher order modulation 
schemes such as 64-QAM are known to be affected particu¬ 
larly worse by PHN. The simulated channel, also employed in 


Table II: summary of the Rayleigh channel parameters 


Carrier Frequency, f c 

5 GHz 

Sampling rate. f s = A- 

25 MHz 

rms delay spread, Tr ™ s 

120ns 

Speed, v 

7 km/h 

Channel taps, L 

to 


1 13H . is given for a normalized amplitude by 1 2711 


H m {k) = Y, c/ ) (r i )e j ^ +2 ^ Dm+2 ^ (31) 

i—0 

where a.i, ipi, fP and r» are respectively the amplitude, phase, 
Doppler frequency and time delay of the z-th propagation path. 
The time correlation for the given channel model is defined 
as 

R t {m -n) = E[H m (k)H*(k)\. (32) 


Assuming uniformly distributed angles of arrival at the mobile 
station, the correlation follows the well known Jakes model 
I 2 H 1 such that R t (m — n) = Jo(2Trf® ax (m — n)), where 
Jo(-) denotes the zeroth order Bessel function of the first 
kind, and f^ lax = ( vf c /c)T is the maximal Doppler frequency 
normalized by the OFDM symbol duration, T. Assuming an 
exponentially decaying power delay profile with normalized 
RMS delay spread r rms , |[ 29ll , the correlation between different 
frequency bins is given by 


R f {k-l) = E[H m {k)H* m {l)] 

1 - e~ L PP 


3 +j2Tr(k-l)) 


(33) 


(1 _ e i/r rms )(i + j2n(k - l)T rms ) 


which satisfies E[\H m (k)\ 2 ] = 1. 

1) Channel fully known at the Receiver: Assuming that the 
frequency response of the channel is known at the receiver, 
the performance of the proposed algorithm is presented and 
compared with previously proposed techniques in 12] and |13>]. 

InH, a Fourier series representation is used to approximate 
the PHN realization at the m-th OFDM symbol, thereby 
aiming to suppress all the influence of the PHN. CPE com¬ 
pensation and symbol estimates are initially obtained using 
pilot subcarriers, after which a DD-FB loop is employed 
for MMSE estimation of the vector of Fourier coefficients 
of the PHN realization which is subsequently used for ICI 
compensation. Similarly, in 11311 , a two stage cancellation 
technique is employed in which the CPE is corrected, followed 
by ICI cancellation using decision feedback symbols. The ICI 
cancellation technique in 130 is derived based on power series 
expansion of the PHN over an OFDM symbol. 

The results in our approach were obtained using Algorithm 
0 but where step 8 and 12 are omitted since the channel is 
already known. The results are shown by the solid curves in 
Figure [8] and Figure [3 for 16-QAM and 64-QAM respectively. 

First, in Figure [8] (solid line), it can be seen that the 
performance of the proposed algorithm without employing 
DD-FB is 0.5dB from the ideal case (no PHN) for a BER of 
10 -4 . On the other hand, for 64-QAM, Figure [9] shows that the 
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performance of the proposed algorithm is 2dB from the ideal 
plot for a BER of 10~ 4 . In both Figures, the proposed method 
shows performance gain in BER compared to the previously 
proposed methods in |H], which employs a single iteration on 
the DD-FB, and the method proposed in UbIi . which is shown 
for 2 and 3 DD-FB iterations respectively for 16-QAM and 
64-QAM. It should be noted that no DD-FB loop is performed 
in the proposed approach. 

2) Unknown Channel Response: For the case where the 
channel is unknown, joint channel estimation and PHN com¬ 
pensation is employed according to Algorithm [2] Performance 
of the techniques in 01 and 01 is also presented for 
comparison. 

In 14], FFT interpolation is employed using the least 
squares estimates of the joint channel frequency response and 
CPE at pilot positions to obtain estimates over all subcarriers. 
The method also proposes estimating the CPE by averaging 
the phase displacements on pilot subcarriers which is used 
to remove its influence from the estimated channel frequency 
response. Previously decoded symbols are then employed 
for least squares (FS) estimation of the vector of Fourier 
coefficients of the PHN realization for ICI compensation using 
previous channel and symbol estimates. Three sets of results 
are presented for the method in 11411 . The first two results 
consider estimation of only one and three Fourier coefficients 
of the PHN realization closest to the DC for ICI compensation 
thus requiring less computational complexity. The third result 
requires inverting a 64 x 64 Toeplitz matrix to estimate the 
entire vector of Fourier coefficients of the PHN realization for 
ICI compensation. 

The proposed method as well as the method in |[13|] use 
D = 3 past OFDM symbols as described in (l24l ) to enhance 
performance of the channel estimation. Performance plots for 
the unknown channel case are shown with the dashed curve in 
Figure [8]and Figure [H for 16-QAM and 64-QAM respectively. 
In both Figure[8]and Figure[9] the performance of the proposed 
technique shows improved performance, in BER, over the 
previously proposed methods in E3, iil. 

In Figure [H the proposed method is 1.5dB from the ideal 
case, (no PHN and known channel response) at a BER of 10 -4 
employing only 2 iterations on the DD-FB loop. On the other 
hand, the method in 11311 . shown for 8 iterations, is 2.5dB from 
the ideal case while the method in 114] is far from the ideal 
case. 

For the 64-QAM system. Figure [9] shows that the proposed 
method provides improved performance compared to ill311 and 
the method in d. With 10 iterations, the proposed methods is 
3dB from the ideal case at a BER of 10~ 3 , while the method 
in El is 5dB from the ideal case for the same number of 
iterations. 

3) Performance with respect to size of K: In Figure 1101 
simulation results are presented showing the BER performance 
of the proposed algorithm for various codebook size K and 
for f)T = 0.01 and j3T = 0.05. The presented result does 
not employ any DD-FB loop. A 16-QAM OFDM system with 
the same description as in Sec. IV-BI is considered assuming 
known and unknown Rayleigh fading channel. It can be seen 
that the BER improves a K grows. However, the relative 
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G3 partial matrix 3 
QJ. full matrix 

■ mi . without DD-FB 
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ICI-cancellation Method in Q] (known channel) 

■Q. 2 DD-FB iterations (known channel) 
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Full channel knowledge and no PHN 
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Figure 8: BER performance of the proposed algorithm compared to B and 
11141] for 16-QAM OFDM. Performance plot for the ideal case (i.e. perfect 
channel knowledge and no PHN) is also given. ICI cancellation technique in 
a is also shown for perfectly known channel. BT = 0.01 



Figure 9: BER performance of proposed algorithm compared to m and B 
for 64-QAM OFDM. Performance plot for the ideal case (i.e. perfect channel 
knowledge and no PN) is also given. ICI-cancellation technique in a is also 
shown for perfectly known channel. BT = 0.005 


improvement as K grows beyond K = 27 is very small. 
For /3T = 0.01, BER improves only slightly at high SNR 
in both cases where the channel is assumed known and where 
the channel is estimated as I\ increases from 27. The relative 
improvement in BER as K increases is more noticeable for 
f3T = 0.05, even though it is still small. 

4) Sensitivity to input parameters: The plots in Figure 
[TT1 show the sensitivity of the proposed algorithm and the 
algorithm in 11311 to the estimated channel and PHN parameters 
relative to the true parameters. The rate at which the PHN 
variance grows is set to BT = 0.03 and the normalized 
maximal Doppler frequency is set to f^ lax = 0.03. All other 
Rayleigh channel parameters are kept the same as in Table QI] 
Moreover, the size of the codebook was K = 27 with Q = 3 
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Figure 10: BER performance of the proposed algorithm on 16-QAM coded 
OFDM for various number of codebook size. 



method in (H (red line) as a function of the input doppler frequency f^ ax 
for various input values of /3T in a 16-QAM OFDM on a Rayleigh channel 
where the true f° ax = 0.03, 0T = 0.03 and SNR = 20dB. 


and J = 4 and SNR = 20dB. Both the proposed method and 
the method in 0 use the past two symbols together with the 
current symbol D = 3 for channel estimation. 

It can be seen that the sensitivity of the proposed algorithm 
as well as the algorithm in 0' to the input values of the 
normalized maximal Doppler frequency, f£ ax is low except 
for f^ ax = 0. However, for higher value of D, the BER 
performance is expected to be very sensitive to the input value 
of frnax- However, it can be seen that the performance of both 
the proposed algorithm and 0 is noticeably sensitive to the 
input values of / 3T. On the other hand, the method in 0 does 
not depend at all on the channel and PHN statistics and thus 
the performance remains unaffected by the estimated input 
parameters. 


5) PLL Oscillator : The above results have been provided 
assuming a free running oscillator in which the PHN is 
modelled by the Weiner process and based on which the 
codebook is derived. In order to demonstrate the applicability 
for a PHN model other than the Weiner process, we employ the 
proposed method using the Ornstein-Steinbeck process [ 30] to 
model the PHN. The Ornstein-Uhlenbeck process is used to 
model the PHN at the output of a voltage controlled oscillator 
(VCO) in a phase-locked loop (PLL) (it]. 

In Ligure [12] a PHN realization, 9(n), for a free running 
oscillator and a PLL is shown. It can be seen that the PHN 
from PLL (Ornstein-Uhlenbeck process) is stable around the 
mean while the PHN from a free running oscillator grows over 
time. 

A detailed analysis of PLL circuit parameters and the asso¬ 
ciated PHN process generated is given in Bill . In Ligure IT2l 
the PHN variance of the Weiner process (free running VCO) 
grows at a rate of (/ 3T) vco = 0.05 while more parameters are 
required to model the Ornstein-Uhlenbeck process. A Charge 
pump PLL is considered which is set to have a 3dB bandwidth 
°f ftp = 20kHz for the low pass filter, f p d = 20kHz for the 
phase detector and f p u = 100kHz for the PLL itself. The PHN 
variance of the VCO and the reference oscillator grow at a rate 
of f3T vco = 0.01 and ( (3T) re f = 2 -9 respectively. Moreover, 
the center frequency of the VCO is f c = 5GHz while the 
center frequency of the reference oscillator is f re f = 100MHz 
and the sampling frequency is f s = 25MHz. A 16-QAM coded 
OLDM with the same description as in Sec. lV-Bl is considered 
assuming unknown Rayleigh fading channel so that channel 
estimation is employed along PHN compensation. 

In Ligure [13] the performances of the proposed method and 
the method in 0 is compared for various input of the PHN 
variance of the VCO, which is denoted by /3T vco , for a free 
running VCO and for PLL at 20dB SNR. Both methods are 
derived based on the Weiner PHN assumptions. The proposed 
method is run in 2 DD-LB loops while the method in 0 
is run in 5 DD-LB loops. The improved in performance for 
PLL (dashed line) compared to free running VCO (solid line) 
can be seen in Figure [13] for both the proposed as well as 
the method in [13]. It can be seen that the lowest BER for 
the proposed method in a PLL is achieved around the input 
/ 3T vco = 0.02 which then starts to increase as (3T vco > 0.02. 
The increasing trend in BER can also be seen in free running 
VCO when /3 T vco > 0.04 for the proposed method. The BER 
performance of the method in [13], on the other hand, seems 
to be stable around 5xl0~ 3 ina PLL and around 9xl0 _3 ina 
free running VCO for an input f3T vco equal to 0.01 and above. 
For the range of input [3T vco values displayed, the proposed 
method outperforms |L3|] in both cases when a PLL and a free 
running VCO is used. 


VI. Computational analysis 

The proposed scheme does incur some computations which 
are not part of the reference schemes 00. These include 
the multiplication with the matrix A ; ) 1 in (f20b . and all the 
other calculations inside the loop over k in Algorithm [2] 
It should be noted that the channel estimation in (l26l >. ( 1281 ) 
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Figure 12: PHN samples, 6(n), simulated as an out from a free running 
oscillator (Weiner process) and a charge Pump PLL (Ornstein-Uhlenbeck 
process). 



Figure 13: Comparison of BER performance of the proposed algorithm and 
the method in fl3ll in a free running VCO and PLL as a function of the input 
(3T vco in which the true (/3T vco ) = 0.05. A 16-QAM coded OFDM and a 
Rayleigh channel where SNR = 20dB used. 


and ( 129) only requires the multiplication of the LS channel 
estimate with a pre-computed matrix. The multiplication with 
A a . 1 is an operation in the frequency domain. More efficient 
would be to do the codebook pre-compensation in the time 
domain and do a separate FFT for each of the K possible 
received signals. In this case, the number of complex multi¬ 
plications to obtain Y h m for k = \..... K is on the order of 
N( 1 + y log 2 N), rather than N(KN + \ log 2 N). A sum¬ 
mary of the number of complex additions and multiplications 
required by the most important equations of our algorithm is 
given in Table [III] 

The method in Cl involves some computations which are 
not needed in our algorithm for instance the successive least 
squares estimation of the coefficients of the power series which 
approximates the PHN in addition to an N x N Toeplitz 


Table III: Summary of total Number of complex operations required per 
OFDM symbol and DD-FB iteration excluding the Viterbi and a few other 
computations. 


Equation 

Additions 

Multiplications 

Complex inversions 

HD) 

KN(N - 1) 

KN 2 

- 

{27) 

- 

KN 

KN 

{26) 

KN(N - 1) 

KN 2 

- 

{25) 

(K + 1) (TV — 1) 

K(2N + 1) + N 

KN +1 

{H 

— 

2KN 

K(N + 1) 

{24} 

KN{N - 1) 

KN 

- 


matrix inversion. The overall number of complex additions and 
multiplications for channel estimation and PHN compensation 
required for one OFDM symbol by the proposed method as 
well as the method in Cl is given in Table IIVI excluding the 
Viterbi algorithm and without taking into account the N x N 
Toeplitz matrix inversion for 131 . The parameter A in Table 
HV1 refers to the order of the polynomial approximating the 


ICI coefficients as presented in 11311 . It can be seen that the 
complexity of both the proposed algorithm and the method 
in 11311 are linearly dependent on K and A respectively. For 
K = 27, A = 6 and N = 64 the proposed algorithm has a 
computational advantage over 113j], This is illustrated in Figure 
[T4l which compares the measured execution time for various 
number of DD-FB iterations between the proposed algorithm 
(Algorithm [2]) and the algorithms in 113] for 64-QAM OFDM 
at 24dB SNR. No DD-FB loop or 0 iteration for the method 
in Cl implies channel estimation with CPE compensation 
which reduces the number of computations significantly. The 
N x N Toeplitz matrix inversion in 11311 has been factored 
out of the measured execution time (since we did not take 
into account the Toeplitz structure). Our simulator is a Linux 
system (Ubuntu 11.04) running on a CPU with Intel-Core 
i7-2600 and uses the IT++ library in C++ for executing 
the algorithms. The results in Figure [14] show that, for the 
chosen parameter I\ = 27, the proposed method provides a 
significantly better performance-complexity trade-off than the 
method in |H. 

The DD-FB stage in H and 112] are not easily parallelized, 
while the most computationally demanding operations in the 
proposed algorithm can obviously be split into up to I\ 
parallel processors or dedicated hardware. The algorithms in 
B5h and 114] also involve an inverstion of an N x N (Toeplitz) 
matrix. However, the method in [14] also presents a solution 
in which only the ICI coefficients close to the carrier need 
to be estimated. This requires solving a less complex linear 
system, e.g. based on tridiagonal matrix algorithm which 
only requires O(N) complex operations. This makes Cl 
less computationally complex to implement than the proposed 
method and m. In addition [14] has the advantage provided 
by avoiding the requirement of having to know the channel 
and PHN statistics. These advantages come at the cost of 
performance in BER as shown in Figure [8] and Figure [9] 


VII. Conclusion 

Based on a simple codebook table which approximates the 
phase-noise statistics by a finite number of realizations, a novel 
PHN compensation approach is introduced. The general idea 
can be applied to a wide range of PHN scenarios. Herein, 
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Table IV: Summary of total Number of complex operations required per OFDM symbol excluding the Viterbi and a few other computations. 


Method No. of Additions No. of Multiplications 

Proposed method (i + 1)(KN(SN — 2) — k + N — 1) (i + 1)(K(2N' Z + 6N + 1) + N) 

[IT] i(AN(N - 1) + N(N - 2)(A - 2)(A - l)/2) i(A(4N 2 + 2N + 1) + A(A - 1)N(2N - l)/2) 

+(i + l)2N(N-1)+(i + l)(2N(N + l)) 



Figure 14: BER performance of Algorithm [2] and the algorithms in 0 on 
64-QAM OFDM as a function of averaged execution time per frame (in sec) 
for K = 27 and A = 6 at 24dB SNR using various number of iterations 
(denoted by i*) in the DD-FB loop. 


we have concentrated on ICI suppression in a convolutional 
encoded OFDM system with scattered pilots. The complexity 
of the algorithm is determined by the size of the codebook. 

Using a codebook of moderate size, K = 27, for a 16- 
QAM uncoded OFDM system in which an AWGN channel 
and a PHN with BT = 0.01, which is was considered, the 
proposed algorithm is 2dB from the ideal PHN free case at 
a BER of 10" 3 . Considering an identical scenario in a 64- 
QAM uncoded OFDM, the proposed algorithm is 2.5dB from 
the ideal PHN free case at a BER of 10 -2 . By comparison, 
the method in Jl3] is 9dB from the ideal case. 

For a coded system, the codebook technique is used along 
with DD-FB loop for combined channel estimation and PHN 
compensation. For a 16-QAM OFDM system over a Rayleigh 
fading channel with BT = 0.01, the proposed method is 
1.5dB from the ideal case (PHN free and known channel) at 
a BER of 1 0” 1 employing 2 iterations on the DD-FB loop; 
while the method in 01311 is 2.5dB from the ideal case using 
8 iterations. For a 64-QAM coded system and PHN with 
BT = 0.005, the proposed method provides a ldB gain at a 
BER 10~ 3 employing 3 iterations in the DD-FB, compared 
to the method in [13] which uses 10 iterations. The gain 
obtained by the proposed method is significant compared to 
the method in [14], which uses a pilot interpolation technique 
for channel estimation. Additional results are also presented 
for both 16 and 64 QAM coded systems assuming known 
Rayleigh channel frequency response at the receiver. The 
proposed method is shown to outperform to the reference 


schemes without having to employ the DD-FB loop. 

The performance of the proposed method with respect 
to the codebook size is also evaluated showing the relative 
improvement in BER as K increases. However, for the given 
example, the relative improvement is very small as I\ increases 
above 27. Simulation results are also shown employing PHN 
from a charge pump PLL, modelled as Ornstein-Uhlenbeck 
process, is also presented with the proposed method showing 
improved performance compared to the methods in El and 
also demonstrating the applicability of the proposed algorithm 
for a PHN process other than Weiner process. 

Therefore, with a moderately small codebook size and lim¬ 
ited number of iterations on DD-FB, an improved performance 
as well as faster execution time is achieved by the proposed 
algorithm. 


Appendix A 

Proof of Equation (IT8] > 

Let us then first consider the case when the number of 
quantization regions, Q , is infinite. In such a case, each 
segment of the selected codebook entry will follow the mean of 
samples of the corresponding segment of the PHN realization. 
Thus the MSE for infinite Q is given by 

J -1 L—l 

MSE (q=oo) = £ £ E e(n) Mn + jL ) - 0(n + jL)) 2 } 

• 7=0 n=0 

J- 1L—l 

= EE E[(0(n + jL) 2 

j =0 n=0 

- 2 8(n + jL)6(n + jL) + 6(n + jL) 2 ] 

= EE E[(9(n + jL) 2 

j= 0 n=0 

-iEEE E[6(n+jL)0(l+jL)] 

j= 0 71=0 1=0 

(N-J)(N + J) _ 2 

- - & e 

6 J 

(34) 

in which the last equality follows from the covariance of a 
random walk process which is E[6(n + £)9(n)] = na 2 for 
l > 0 . 

When Q is finite, an additional error is introduced due 
to the quantization of segment averages. As in the design 
of the codebook, we assume that -0 is set to the average 
of the first segment and thus no quantization error occurs 
in the first segment and the PHN realization crosses the 
first sample of each segment of the trajectory. Additionally, 
by assuming that this error in each segment is identically 
distributed and independent of the difference between the PHN 
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and the codebook trajectory, the MSE due to accumulated 
quantization error in each segment is La 2 , where L = N/J 
and a q is the standard deviation of quantization error. The 
MSE due to the quantization error over the entire trajectory 
is therefore, MSE(q <00 ) = (J — l)La 2 , where (J — 1) is due 
to the assumption that the first segment has no quantization 
error. The approximate MSE of the selected trajectory as an 
estimator of the PHN realization then becomes 


MSE = {N + ^ J) a 2 t + L(J - l)u g 2 . (35) 


In order to determine a 2 , we assume that the PHN realization 
is exactly on one of the codebook entries in the last sample 
of the previous segment. In this case, the quantization error of 
the j-th segment average is equal to the quantization error of 
the j-th increment, i.e., no influence of previous increments. 
The variance, n 2 , of the quantization error is then given by 


Q 

a l = Y, E ^ X - f ) 2 ! X e Ki\P{X e Hi)) 

2 — 1 
Q 

= V(X G Hi){E[X 2 \X G Hi] - 2xiE[X\X elZi]+x 
2—1 
Q 

= J2 V ( X e Ki)(E[X 2 \X G Hi] - X?) 

2=1 

(36) 


where the last equality follows from E[X\X G 7?.,] = 
Xi which is given by (fl6l >. The conditional Expectation 
E[X 2 \X G Hi] is given by 


E 


X 2 |A' g Ri 


J R . x 2 f x (x)dx 
V[X G 


(37) 


where P(I G K,) = f R fx(x)dx. When fx{x) a Gaussian 
function 


J x 2 f x (x)dx = al/2^evf(x i+1 /V2a x ) - er f(xi/V2a x ) 

— a x V2jt^Xi+\ exp(—0.5 x 2 i+l /a 2 x ) — Xi exp(— 0.5x 2 /a 2 ) 


and 


P(XeH i )= [ fx(x)dx 

jRi 

= 0.5^erf(x i+ i /V2a x ) - erf(x;/-\/2cr x )J . 


(39) 


Equations (l38l) . l39l and (IT6l > are then used on (l36l > to obtain 
the variance, a 2 , of the quantization error. 
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